Pulsed Laser Cooling for Cavity-Optomechanical Resonators 
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A pulsed cooling scheme for optomechanical systems is presented that is capable of cooling at much faster 
rates, shorter overall cooling times, and for a wider set of experimental scenarios than is possible by conventional 
methods. The proposed scheme can be implemented for both strongly and weakly coupled optomechanical sys- 
tems in both weakly and highly dissipative cavities. We study analytically its underlying working mechanism, 
which is based on interferometric control of optomechanical interactions, and we demonstrate its efficiency 
with pulse sequences that are obtained by using methods from optimal control. The short time in which our 
scheme approaches the optomechanical ground state allows for a significant relaxation of current experimen- 
tal constraints. Finally, the framework presented here can be used to create a rich variety of optomechanical 
interactions and hence offers a novel, readily available toolbox for fast optomechanical quantum control. 



Introduction - Micro- and nanomechanical resonators are 
currently emerging as new quantum systems f]^. Their inte- 
grability in a solid state architecture offers attractive opportu- 
nities for quantum information objectives such as mechanical 
quantum registers J^lsl. optomechanical quantum transduc- 
ers (01 or quantum memories Js]]. At the same time, their 
size and mass promise access to a hitherto untested regime 
of macroscopic quantum physics jg-lsll. The field of cav- 
ity quantum optomechanics |T3-13D utilizes methods from 
quantum optics in combination with optomechanical radia- 
tion pressure interactions to achieve this and experiments are 
progressing r apidly - the strong coupling regime has been 
demonstrated lll9l - l2lll . and optomechanical analogues of elec- 
tromagnetically induced transparency n22[ 12311 have demon- 
strated first steps towards mechanical storage of light. 

A prerequisite to achieve full coherent control over me- 
chanical quantum states is to operate these systems close to 
their quantum ground state and to achieve coupling rates that 
exceed all other decoherence rates. Optomechanical cooling 
close to lfl4Ul6lEllE3ll and even well into the quan- 

tum ground state of micro- and nano-mechanical devices have 
been realized. However, most mechanical devices are intrin- 
sically connected to a hot environment through their supports, 
which results in large heating rates. Efficient cooling therefore 
requires to minimize the thermal coupling, either by operating 
in a cryogenic environment lll4l - ll8l 12411 or by decoupling the 
mechanical resonator from its environment 11251 - 127 1 . 



As a consequence, a cooling scheme that can beat the me- 
chanical heating rate for systems that are not cryogenically 
cooled or mechanically decoupled from the environment is 
highly desirable. For the most widely used scheme, sideband 
cooling 128-311, the cooling rate F, and conversely the time 
required to approach the ground state, are inherently limited 
by the mechanical frequency, v. Specifically F < v, due to use 
of the rotating wave approximation (RWA). Recently it was 
demonstrated in the context of ion trap physics that pulsed 
schemes can break the speed limit set by the oscillator fre- 
quency II32II . by using interference between optical pulses in- 
cident on the system that is being cooled. This generates an ef- 



fective cooling (red-sideband) term with large couplings, thus 
avoiding the RWA limitations. Unfortunately, this approach 
cannot be directly used as the nature of the cavity coupling 
is significantly different (generally non-linear for the optome- 
chanical system, and when linearized, h.o. (harmonic oscil- 
lator) to h.o. coupling v. h.o. to spin coupling for trapped 
ions) as is the nature of the subsystem from which energy is 
eventually removed (optical cavity with a wide range of qual- 
ity factors, v.s. the ion with a finite Hilbert space and fixed 
physical properties). Moreover, one needs to overcome the 
instability issue in optomechanical systems, which does not 
exist for trapped ions. 

In this letter we demonstrate how a sequence of fast pulses 
adds a term to the effective optomechanical interaction Hamil- 
tonian which approximates the cooling (also known as beam- 
splitter, anti-Stokes and red sideband) operator x,„Xc + PmPc °^ 
ab'' + a^b. Here a is the annihilation operator of the cavity and 
XcPc are its quadrature operators, b is the annihilation oper- 
ator of the mechanical oscillator with corresponding quadra- 
tures x„,,p„,. The technique is shown to be experimentally fea- 
sible in both the good cavity (k <k v) limit, where it is capable 
of reaching the ground-state much faster than the oscillator 
frequency, and the bad cavity limit (k > v), where sideband 
cooling is incapable of approaching the ground state. 

Physical Setting - Cavity optomechanical systems are mod- 
eled as an optical cavity field which couples to a mechanical 
resonator by way of radiation pressure. The optical (mechan- 
ical) mode, oscillating at a frequency w (v) is characterized 
by a relaxation rate k (ym). The optical mode is driven by a 
detuned laser field of frequency w; with a strength Q. The 
corresponding Hamiltonian is 

HoM = Aa^'a + vb^'b+^a^'a [b' + Z7)h-Q (a'c""* + ae'''') (1) 

where (/> is the initial driving phase, A = w-w; is the cavity de- 
tuning w.rt the laser frequency, and go is the optomechanical 
coupling rate. 

Linear Approach - For continuous driving, cooling is 
achieved by invoking the RWA. In order to derive a pulse se- 
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quence which generates cooling directly via the cooling oper- 
ator introduced above, i.e. without RWA, we will need to an- 
alyze the dynamics of the system. Assuming, for now, weak 
coupling, go <*: v, one may avoid the complexity of the non- 
linear nature of the interaction by considering a linearized ap- 
proximation of the system. 

As this scheme makes use of rapidly changing driving of 
the cavit y m ode, the usual Hamiltonian linearization proce- 
dure 

HI 12811 cannot be trivially applied to eq. ([T]i. Rather, we 
move to a frame co-moving with the state of the cavity, by 
applying a time-dependent canonical transformation (see Ap- 
pendix A for additional details). In this frame, a is redefined 
as a small perturbation, allowing us to replace the quadratic 
coupling with a linear one. This results in the Hamiltonian: 



Hii„, = Afl"' "fl + vb^'b + (G (t) a + G* it) fl"'") x„, + \G (f)l^ 



(2) 



with G (f) = igoe-"-^-'"'" Q (f) e'^^-''^'' dt' . 

The linear nature of the dynamics now allows us to re- 
phrase the dynamics as an equation-of-motion of the covari- 
ance matrix ll34ll . We define a vector of quadrature opera- 
tors R = (xc, Pc, x„t, p„,y and the covariance matrix as jij = 
2Re {{RiRj^ - {Ri) (Rjfj- The corresponding equation of mo- 
tion is jj - My + y^M^^H-lfwithM - SV-^Pwhere Vis 
the potential matrix (///,„. - R^VR), S the symplectic matrix 
andP = diagil, 1,0,0). 

Denoting a series of « control pulses (i.e. n modulations 
of the cavity driving laser), {Hc\, . . ■,Hc„] of corresponding 
durations (fi...f„) (with free evolution implicit), and using 
the Baker-Campbell-Hausdorff (BCH) jssll equivalency, one 
may compute the equivalent control Hamiltonian He from 
gTt„{Ho+H„,)_ ,^^i,(Ho-i-H,.,) _ 'i)(Ho+w,.)^ where Hq as in 

eq. (H for G = 0. 

The equivalent control Hamiltonian to the control sequence 
{-GxcX„,, +GpcX,„, -GpcXm, -GxcXin] for (f2, f i , f i , fa), with 
A = V (to second order in the BCH series) is 



Hc=G (ff V - 2t2) XcX,„ - Gt\vpcPm- 



(3) 



Setting f2 = ^iV, the desired -Gt\v{XcXm + PcPm) cooling op- 
erator is achieved. While higher-order elements in the BCH 
series contribute terms undesirable for cooling (e.g. squeez- 
ing), both analytical estimates and numerical studies show 
these disruptions to be well-contained for pulse durations 
shorter than the mechanical frequency. Moreover, as will 
be shown later, optimal control can further reduce undesired 
terms, improving upon sideband cooling w.rt both the achiev- 
able final temperature and the required time. Finally, it is im- 
portant to stress that the BCH operators are a two-edged sword 
- they create the cooling operators from commutation relations 
between the driving terms and the free Hamiltonian, but at the 
same time they transform the dissipative elements (Lindblad 
terms) in a similar fashion. As a result, stronger driving en- 
hances both cooling and dissipation terms (modified by the 
BCH relations), resulting in sub-linear advantage of such an 
approach to further increase cooling rates. 



Optimal Control - We have examined the performance and 
limits of the proposed cooling sequence using search methods 
from optimal control. We took multiple approaches, includ- 
ing (a) two stage optimization - initially optimizing the pulse 
amplitude and subsequently optimizing both amplitude and 
phase; (b) using the analytically derived sequence in eq. (O 
as an initial point of the optimization; (c) in case of strongly 
dissipating cavities, a "telescoping" series of optimizations is 
used, slowly increasing dissipation, with the sequence result- 
ing from optimization k serving as the initial condition for 
optimization ^ + 1 (d) a similar series of optimization is used 
to gradually shorten overall cooling times (e) random starting 
conditions and simultaneous optimization of all control pa- 
rameters. Fig. ([TJ shows the results of multiple optimizations 
for varying values of k, and a comparison to the results of the 
analytical sequence before numerical optimization. Fig. ^ 
shows optimization results for the bad cavity regime. There, 
we allow the system state to become squeezed alongside the 
desired cooling and optimize for the squeezed phonon count 
(measured by the reciprocal purity iIstIIssI '). 

An example of such an optimal sequence has been obtained 
using the full range of interactions, (7?e[G]jCc + Im[G]pc) x,„. 
Starting with an initial thermal phonon occupation of 100, and 
for K - 0, we are able to achieve a final occupation below 
2 X 10"^ in less than 0.62;7rv ' (see Appendix D). Another ex- 
ample, illustrated in figure ([3) presents the detailed behavior 
of a highly dissipative cavity k = 2167v. For such settings, it 
is much harder to suppress all undesirable terms in the higher 
orders of the BCH series. However, the system can be cooled 
in under 10"^^ to less than one (squeezed) phonon, if, again, 
we allow the system state to become squeezed. All optimiza- 
tions have been performed using QLib 135 1 . 

Non-Linear Case - For systems with large coupUng go 



1 39l |40tl it is necessary to treat the full nonlinear interaction as 
presented in eq. ([TJ. Unfortunately, in the case of the single- 
cavity system analyzed thus far, it does not appear to be pos- 
sible to generate a cooling operator using a simple sequence 
of pulses analogous to that used for the linearized system (see 
Appendix B). 

To overcome this limitation, a double cavity device, where 
the oscillator is coupled to two identical cavities (on the left 
and right of the oscillator), is proposed ll4ll - l43ll . This geomet- 
rical arrangement results in the coupling of the two cavities to 
the oscillator to be of opposite sign, which allows one to can- 
cel the mechanical oscillator driving using opposing radiation 
pressure of the two cavities. The undriven system Hamilto- 



Ho - A(a\ai + a\ai) + vb^b + ^^{a\a\ - alfl2)(fe' + b), (4) 



V2^ 



where the subscript distinguishes the two cavities and A is 
the detuning of both cavities. We define the symmetric and 
antisymmetric modes, aa - -^(ai + aa), a.? = "^("^i ~ '^2) 
and the dimensionless quadrature operators x = -^(a' +a^, 
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FIG. 1. Results of several covariance matrix optimization proce- 
dures. All with an initial phonon thermal occupation of 10, |G| < lOv, 
7 = and taking at most 0.8 x ^. The red and magenta data sets in- 
volve optimizations with random initial pulse-sequences, optimized 
with partial (xrJi:,„-only) and full coupling respectively; the random 
initial settings and nature of global optimization is responsible for 
the poor performance at /<• « 0.03. Final temperatures achieved by 
sequences of 30 pulses (based on the four pulse analytical formula- 
tions, repeated 7.5 times) are shown before and after optimization 
(blue and cyan). Long cyclic sequences (75 repetitions of the 4 pulse 
analytical sequence) are shown before and after optimization and are 
represented by the dark and light green sets. Application of the last 
set, pre-optimization, for k = 0, appears at the left axis of the plot. 
Note that as often with numeric optimization, the points may repre- 
sent local optima. 
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FIG. 2. Performance of optimized cooling sequences for k > v. 
Starting with the analytically derived sequences, a series of optimiza- 
tions was used to generate cooling sequences for successively higher 
dissipative systems, with the result of one optimization serving as the 
starting point for the next. Note that in this case, the resulting system 
state is squeezed, as a result of higher-order BCH terms. 



p = - aj, re-expressing the Hamiltonian as 

Hq = A(fl>a + ajflj) + Vb'b + go(XaXs + PaPs)Xm- (5) 

Building the proposed cooHng sequence in stages, let us 
first examine the terms generated by a {-Q.Xa,Q,Q.Xa} se- 
quence, goQ.tipsX,„ + 2ls.Q.t\pa. By driving the cavity during 
the free evolution period with a pa pulse, the coefficient for 
Pa, above, can be controlled at will. By defining this addi- 
tional pulse as -{a + 2)AQ.tipa the generated terms transform 
to: 

HnLI = gO^tiPsXm - aAiltipa. (6) 

The coohng sequence for the double-cavity setup is 



{{-Q.Xa,/3pa, QXfl), 0, {Q.Xa,/3pa, 



(7) 



where the nested notation emphasizes the nesting of se- 
quences used. Here /3 = -{a + 2)AQfi and the pulse dura- 
tions are ((fi, f/, fi), f^, (fi, f/, fi)) with t2 = ti + tf + t\ and 

Assuming Q.ti » 1, the free evolution implicit in Hmi of 
eq. ^ can be neglected; the inverted sequence {QjCq, 0, -Q.Xa] 
together with the additional driving during the free evolution 
will yield eq. (|6]l with an inverted sign for Q., resulting in what 
can be viewed as a sequence-of-sequences, i.e. a nested pulse 
sequence. Defining ti = Qfif2, we get the terms generated by 
eq. O to be 



Hnli = +(« - 2)goATiXj.x„, -I- 2govTip,„ps 



+'^8qTiPs - 3AgoTiX„ - TigoXaX,„ 

+2aTiA'Xa. 



(8) 



When A - 2v/(a - 2) this Hamiltonian contains the cool- 
ing operator for the mechanical oscillator using the symmet- 
ric mode of the cavities, i.e. the XsX,„ and psPm terms appear 
with identical coefficients. Setting a - 4- one may meet the 
resonance condition A = v, which is important in some exper- 
imental layouts, as stronger interactions are achieved at reso- 
nance. An additional Xg driving pulse, of area -2aTitiA^, can 
be added to the sequence to counter the last term in eq. (jSjl, 
provided the correction term and the entirety of the sequence 
in eq. ^ are combined via the Trotter decomposition 113 611 . to 
suppress any high-order BCH terms. The two quadratic terms 
can be absorbed into the corresponding frequency terms by 
means of a Bogoliubov transformation. 

To achieve ground-state cooling in a time shorter than 
the mechanical oscillation, following eq. (O, we require the 
beam-splitter prefactor Q.ga » max (v, /<•)", choosing the f,'s 
to all be one order of magnitude smaller than l/max(v, /c), 
to suppress higher-order BCH terms. Finally, note that when 
G ^ V one may do without the periods of dedicated free evo- 
lution, although pre-factors will be different. 

Experimental feasibility - To obtain both a cooling rate F 
(defined as the pre-f actor of the cooling operator in the Hamil- 
tonian) beyond the limitation of continuous sideband-cooling. 
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FIG. 3. Detailed behavior of highly optimized cooling pulses for an 
optomechanical system with k = 2167v, with 10 phonons in the ini- 
tial state. Top panel contains the reciprocal purity for the mechanical 
oscillator and cavity; second panel the quadrature squeezing; third 
panel the covariance matrix elements and the forth panel details the 
driving. Note this is a regime in which sideband cooling is inappro- 
priate, as the overall cooling time is orders of magnitude below v. 



intrinsic limits of conventional continuously pumped cooling. 
Our scheme is based on generating the cooling interaction, by 
quantum interference of successive pulses. While already a 
simple analytical approach provides otherwise unachievable 
cooling rates F > v, the use of optimal control methods can 
further enhance these rates. We have also shown that cur- 
rent optomechanical configurations could achieve dramatic 
improvements in their experimental performance. Following 
the methodology presented in this work, it is possible to gen- 
erate a rich class of optomechanical interactions, for exam- 
ple the down-conversion (two-mode squeezing) interaction or 
various non-linear terms. This establishes a new and complete 
tool kit for fast preparation and manipulation of optomechani- 
cal quantum states and may very well provide a route towards 
room temperature quantum optomechanics. 

This work was supported by the AXA Research Fund, 
by the EU STREP projects HIP, PICC, MINOS and EU 
project QESSENCE, by the Austrian Science Fund (FOQUS, 
START), by the European Research Council (ERC StG), as 
well as by the Alexander von Humboldt foundation. 



While finishing this paper we learned of 14511 . which also 
treats pulsed cooling schemes for optomechanical systems. 



i.e. F > V, and to reach ground-state cooling faster than the 
mechanical frequency, readily available experimental param- 
eters are sufficient: for an optomechanical Fabry-Perot cavity 
one can easily obtain v - 27rl0*' Hz, nieff = 5 x 10"" kg and 
K - 0.75v ( F5, T^l, which yields go - 15 Hz ^ v and hence 
satisfies the linear regime of pulsed laser cooling (assuming 
a cavity length L = 10"- m and an optical pump wavelength 
A = 1064 nm). 

Previously achievable cooling rates have been of the order 
of F Si 10"' XV lllSn . With our method, and making use of 
optimal control, we fine-tune a 10-pulse cooling sequence of 
total duration 0.75 =f and pulse energies of a: 40 nJ per pulse, 
which can be created directly from amplitude modulating a 
Q.5W continuous-wave laser beam, to obtain a cooling rate 
F = 1.3 V. For an initial temperature of T = IK (rii » 2900 
phonons) and y,,, 300Hz (Q - 2 x 10"*) this sequence will 
reach mechanical thermal occupancies on the order of ^ 

0. 1. 

As a second example we consider a highly dissipative op- 
tomechanical Fabry-Perot double microcavity with y = 2;rl0^ 
Hz, nieff - 10"'" kg, K - 2 X 10^ v and individual cavity 
lengths Li - L2 - 4/1, as has been suggested in 14411 . The 
resulting go - 10'' Hz > v satisfies the nonlinear regime of 
pulsed laser cooling. Our method requires a set of laser pulses 
of length <*: 50 ps with a maximal peak power of 1 kW, which 
is available, for example, in the form of Q-switched lasers. 
This will result in a net cooling rate F = 10'' v. For y„, k 1, 

1. e. a Q-factor of 6 x lO"*, this would already allow cooling to 
the quantum ground state starting from room temperature. 

Conclusions - We have introduced a novel pulsed cool- 
ing method for mechanical oscillators, which surpasses the 



APPENDIX A - THE LINEARIZATION PROCEDURE IN 
THE STRONG-DRIVING LIMIT 

Starting with a simplified view of the cavity (limited to X 
driving only) 



H = A(a''fl) + Q.{t){a + a~) 

we derive the expectation value a of a{ty. 

a{t) - i[H, a] - -iAa - /Q (f) 
a - -iAa - /Q (f) 



(9) 



a{t) - aoe - ie 



r Q.{t')e'^''dt' 
Jo 



(10) 



= aoe-"" -fit) 



We apply a displacement, D - e" """" , to a state-vector 
I^F) = D I*!*), and derive the transformed Hamiltonian H by 
expanding the LHS of the Scrodinger equation for the dis- 
placed vector, ll^P) = -iH\'^). The result is {a* a - aa')\^) - 
iDHD^\'¥} = -iHl'P}, which impHes 



H = DHD ^ + i [a* a - ' ) 

= A (fl"' fl) + \a\- A + 2Q it) {a + a*) 



(11) 



Turning our attention to the interaction term (which we ig- 
nored in the discussion above), one can show the coupling will 
be proportional to 



g(a^ aj + g (aa ' + a*aj 



(12) 
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Finally, note that the exercise above can be repeated for the 
dissipative system, replacing eq. ( fTOl i with a = -Ka - iAa - 
/Q (?) and deriving the modified Liouvillian L. 



APPENDIX B - ANALYTICAL APPROACH TO 
NON-LINEAR SINGLE-CAVITY COOLING 



Let us analyze the eff'ect of applying the pulse sequence ap- 
propriate to the linearized system to the non-linear case. As 
the pcXin pulse is not directly available in the non-linear case, 
it must be generated by a BCH subsequence {-Oa>, Ho, 
Unfortunately, this sequence also introduces a mechanical 

driving term of the form go (fifp) Xm that cannot be countered 
directly by cavity driving, but must be removed by displacing 
the mechanical operators. 

Operators involving both the mechanical oscillator and the 
cavity present another issue, as they generate terms cubic in 
the position quadrature of the oscillator and render the system 
unstable. Specifically, when nesting BCH operations, i.e. ap- 
plying for a time fi the pulse IgoOtpPcX,,, generated by the 
subsequence, in order to obtain the beam-splitter operator. 
The BCH relations imply the substitutions — > Xf - Q.Xm 
and p„j p„ + Slpc- Defining Q' = IgoOtp the effective 
Hamiltonian then reads: 



H ^ Ho + 2vD.'tipcp,„ - 2AD! t]_XcX„ 



+v{Q.'hf pI + A{Q!tiY X, 



2 J2. 



(13) 



2 „3 



-2gQQ.'tiXcxi+ go (Q'fi) x: 



The terms in the first line correspond to the beam-splitter op- 
erator when A = V. The second line presents quadratic terms, 
for which a Bogoliubov transformation would be in order to 
absorb them to the frequency terms of the original Hamilto- 
nian. The cubic terms in the last line are proportional to go- In 
the case of systems with go <s: v, these terms can be neglected 
and hence the linear approximation survives. 



APPENDIX C - CAVITY DYNAMICS 



The detuning A - a> - a>c, with a>c - f c, provides us with 
the only relevant phase difference: 



£,„ (f + t) = Re'^^Ei,, it) + i ^Eou, (f + t) . 



(15) 



Assuming that the driving laser is almost in resonance with 
the cavity eigenfrequency (At «: 1) the phase term can be 
Taylor expanded to first order: 

£,„ (f + t) = (1 + /At) Ei„ (t) + i y/fE,„„ (t + r). (16) 

For lossless cavities R - I - T x 1, and therefore TAt terms 
can be neglected: 



Ei„ (f + t) = Ein (t) - TEi„ (t) 



+ iATEin{t) + i<fEou,{t + T). 



(17) 



In the limit where the dynamics in the cavity are negligi- 
ble during the short r timescale, a diff'erential equation can be 
derived: 



Ein (f + t) - Ein it) T 

= Ei„ (t) 

T T 



(18) 



Since f = 7: 

dEi, 



dt 



+iAEi„ (0 + i Eout {t + T). 

T 



-KEin (t) + iAEin (0 + i yf^—Eo 



(19) 



It is possible to relate this equation to the regular definition 
for the driving: 



(20) 



The power P of an electromagnetic wave can be further ex- 
pressed in terms of the field amplitude E as P - ^"l^^f — 

dt^"^' ~ fo£"^c, where A is the area of the cross section of 
the laser and c the speed of light. Thus: 



(21) 




A differential equation for the amplitude of the field in the 
cavity £„, (f) as a function of a variable external field Egut (0 
of frequency <y is derived here. First of all, a recursion for- 
mula can be found by considering the state of the field inside 
the cavity after a round trip time j = —, with L the length 
of the cavity. It can be considered to be the sum of (a) the 
previous field after bouncing off the mirrors of each end of 
the cavity (of reflectivity coefficient V^) and having picked 
up the phase e"^^ and (b) the field that has entered the cavity 
during this time through the front quarter-wavelength mirror 
of transmissivity Vt. This can be expressed as: 



£,„ (f + t) = Re'^^'Ein (t) + i ^E„„ {t + r). 



(14) 



The pre-factor can be expressed as = so 

that it can be interpreted as the inverse of the electric field 
associated with one photon. The eq.(fT9l) can be rewritten as: 





Ein (t) 



1 fia> 

+ iAJ—Ei„(t) + in 



(22) 



which is by definition: 
da 



— — —Ka + iAa + iQ. 
dt 



(23) 
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0.05 0.1 0.15 0.2 0.25 0.3 035 0.4 0.45 0.5 0.55 



t [2JI/V] 



FIG. 4. Sample pulse sequence optimized for the full linear interac- 
tion, with initial phonon occupation of 100, G,„oi = v, achieving a 
final occupation below 2 x 10"'. 



with a the field amplitude inside the cavity. 

An alternative route to obtain the same result is to transform 
the square of the last term in eq.(fT9]l into; 



EL n it) £,2 



~^ottt _ y-y — QUI 2^ ^ — QUI _2 — QUI (24) 



U ] P 



T n(t) T Ho) n (t) Ho) n (f) 
So the equation has to be rewritten as: 



—— = -KEi„ + lAEin + — 

at \ Tiu y/n 



(25) 



In order to get an equation for the field amplitudes, one has to 
divide by the field amplitude carried by one phonon (^): 



da 
dt 



-Ka + iAa + i ■ 



(26) 



APPENDIX D - SAMPLE OPTIMIZED PULSE SEQUENCE 
See fig. (Ell. 
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